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Visualization and Interpretation of Attosecond Electron Dynamics in 
Laser-Driven Hydrogen Molecular Ion using Bohmian TrajectoriesEl 
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Wc analyze the attosecond electron dynamics in hydrogen molecular ion driven by an external intense laser 
field using ab-initio numerical simulations of the corresponding time-dependent Schrodinger equation and 
Bohmian trajectories. To this end, we employ a one-dimensional model of the molecular ion in which the 
motion of the protons is frozen. The results of the Bohmian trajectory calculations do agree well with 
those of the ab-initio simulations and clearly visualize the electron transfer between the two protons in the 
field. In particular, the Bohmian trajectory calculations confirm the recently predicted attosecond transient 
localization of the electron at one of the protons and the related multiple bunches of the ionization current 
within a half cycle of the laser field. Further analysis based on the quantum trajectories shows that the 
electron dynamics in the molecular ion can be understood via the phase difference accumulated between the 
Coulomb wells at the two protons. Modeling of the dynamics using a simple two-state system leads us to an 
explanation for the sometimes counter-intuitive dynamics of an electron opposing the classical force of the 
electric field on the electron. 

PACS numbers: 33.80.Rv 33.80.Wz 
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I. INTRODUCTION 

The causal interpretation of quantum mechanics by de 
Broglie and Bohm provides the concept of trajectories for 
the dynamics of microscopic objectsi*^. These trajecto- 
ries, called Bohmian trajectories or quantum trajectories, 
are navigated by the wavefunction. Conversely, if we re- 
gard the probability density (i.e., the squared modulus of 
the wavefunction) of the system as a fiuid, the fiow of this 
fiuid can be elucidated by the quantum trajectories^. It 
is this characteristics of the Bohmian trajectories that we 
want to utilize in this article to visualize and analyze a 
recently revealed counter-intuitive electronic motion in 
molecular ion exposed to intense laser light on an 
attosecond time scaled"—. 

Previously, we found in ab-initio numerical simulations 
that at intermediate internuclear distances (i.e., be- 
tween the equilibrium distance and the dissociation limit) 
in an intense infrared laser pulse shows multiple bursts 
of ionization within a half-cycle of the laser field oscilla- 
tion^. This ionization dynamics contradicts the widely 
accepted picture of strong-field ionization, namely that 
an electron leaves the atom or molecule with largest 
probability at the peaks of the oscillating electric field 
of the laser. For example, in the often used tunnel ion- 
ization picture the electron tunnels through the barrier 
created by the binding potential of the ionic core and 
the electric potential of the laser field. This tunnel bar- 
rier is, of course, thinnest when the electric field strength 
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is strongest, which leads to the above mentioned expec- 
tation for the most likely time instant of the electron 
escape. We identified that the multiple ionization bursts 
are induced by a previously reported ultrafast transient 
localization of the electron density at one of the protons^. 
This attosecond dynamics can cause that the electron 
density near the tunnel barrier in the molecular ion is 
highest when the external field strength is below its peak 
strength. Correspondingly, the electron does not tunnel 
most likely at the maxima of the field but at other time 
instants. 

We confirmed and extended earlier interpretations^ 
that the internal dynamics of the electron is a result of 
a strong and exclusive trapping of the population within 
a pair of states of opposite parity, so-called charge reso- 
nant states^. Consequently, the attosecond localization 
dynamics of the electron in the hydrogen molecular ion 
driven by the laser field can be successfully reproduced 
using a two-state modet. The dynamics can be also un- 
derstood by the phase-space fiow of the electron probabil- 
ity density regulated through so-called momentum gates 
which are shifted in time by the vector potential of the 
external laser field^. 

However, results of numerical simulations for the flow 
of the electron probability density often do not reveal 
many details and cannot provide much further insights 
into the internal electron dynamics in the molecular ion. 
We therefore use the concept of Bohmian trajectories to 
provide a complementary picture of the dynamics. By 
analyzing the motion of the trajectories we furthermore 
clarify that it is the phase difference of the wavefunction 
at the two protons that is the origin of the force which is 
sometimes driving the electron in the direction opposite 
to the strong electric field of the laser light. 

The rest of the paper is organized as follows. In SecHIl 
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wc present the model for used for our analysis. In 
Sec. Illli we present the results for the Bohmian trajecto- 
ries and compare them with those for the electron prob- 
ability densities obtained from ab-initio numerical simu- 
lations. We identify the phenomena of transient electron 
localization and multiple ionization bursts in the time 
evolution of the trajectories. In Sec. IIVI we make use of 
the Bohmian trajectory calculations to provide an analy- 
sis of the sometimes counter-intuitive electron dynamics 
in the hydrogen molecular ion. Finally, Section |V] con- 
cludes the paper. 



II. THEORY 

In recent studies of laser induced dynamics of the 
hydrogen molecular ion the full three-dimensional elec- 
tronic motioni^ or the two-dimensional electronic mo- 
tion in cylindrical coordinates along with nuclear motion 
in one dimension have been taken into account^"—. In 
the present study, we consider a simpler onc-dimensional 
(ID) model of in which the positions of the protons 
arc fixed in space. We have shown that the electron lo- 
calization dynamics inside the molecular ion as well as 
the phenomenon of multiple ionization bursts does not 
change in higher dimensional modela^^i^. In particular, 
we found that the nonadiabatic coupling between the 
electronic and nuclear motions is not essential for the 
internal electron dynamics. 



A. ID fixed-nuclei model of 

In our ID model of with fixed positions of the pro- 
tons, the internuclear axis was assumed to be parallel to 
the polarization direction of the linearly polarized laser 
light. The Hamiltonian of this model system is given by 
(Hartree atomic units, e = m = h = 1, are used through- 
out this article unless noted otherwise) 



H{t) = -l-^ + Vciz;R) + Vi^{z,t), (1) 
2 az^ 



where z is the electron position measured from the center- 
of-mass of the two protons. The Coulomb interaction be- 
tween the electron and the two protons was approximated 
by the soft-core potentia l^^'^^ . 



Vc 
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yJ{z + R/2Y+a ^{z-R/2y + a 



(2) 



where a is the soft-core parameter. The laser-electron 
interaction was expressed in the length gauge as 



Vi^{z,t) ^ zE{t). 



(3) 



The laser electric field E{t) is related to the vector po- 
tential A(t) by 



dAjt) 
dfA{t) 



dt 



sin(a;i -|- (/?) -I- fji{t)uj cos(cjt + ip) 



(4) 



where we used the following form of the vector potential: 



Ait) --f Ait) sin 
fA{t) 



Aosin2(^) (0<t<T) 







(otherwise) 



(5) 
(6) 



The full-width at half-maximum (FWHM) of this enve- 
lope function, /^(i), is equal to T/2. 



B. Propagation of the wavefunction and the quantum 
trajectories 

With the Hamiltonian given as above, the wavefunc- 
tion ^'(z, t) was propagated according to the correspond- 
ing time-dependent Schrodinger equation (TDSE), 



i^^^iz.t) = Uit)^iz,t). 



(7) 



This TDSE was solved numerically using the second- 
order split-operator method on the Fourier grid^^rJ^. 
The spatial and temporal grid intervals used for the sim- 
ulations were Az = 0.152 and At = 0.0245 or smaller. 

At the same time as the wavefunction 5'(z,<) was 
propagated in time, the Bohmian trajectories {zj(t)|j = 
1, . . . , iVtraj} were propagated as well by solving the equa- 
tion of motion. 



^ = i;(zj(0,t). 



(8) 



where the velocity field v(z,t) is given by the phase gra- 
dient of the wavefunction. 



*(z,t) = i?(z,t)exp(iS'(z,t)) 
with i? > and 5 G M, as 

dS 



viz,t) 



dz 



(9) 



(10) 



The following identity was utilized in the actual compu- 
tation: 



ds 

IT = ^^^^ 
oz 



1 



(11) 



The ordinary differential equation ([5|) with respect to t 
was solved by the fourth order Runge-Kutta scheme^S 



3 



with the fixed step size of 2 At. i.e., twiee the time step 
of the wavefunction propagation. The wavefunction value 
at every other step of its propagation was used to evaluate 
the velocities of the trajectories at the mid-point of one 
Runge-Kutta step to achieve the fourth order accuracy. 
We may note parenthetically that expression ([T0|) for the 
velocity field is valid for the wavefunction in the length 
gauge representation used in the present study. In the 
velocity gauge representation, the velocity field is given 
by v{z,t) = dS/dz-A[t^. 

The initial positions of the quantum trajectories were 
distributed at a regular interval, 5z, and for each trajec- 
tory we assigned the weight. 



o "55 
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dz|-0(z,to)P, 



Oj{to) 



%(io) 



5z 5z 
zjito) ~ — <z < Zj{to) + — 



(12) 
(13) 




At the limit of 5z — > 0, we may consider that the weight 



FIG. 1. Comparison of the quantum trajectories (a) and the 
electron probability density (b) for the ID model of at 
R = 7 for the interaction of with an electric field (c) of a 
laser pulse (peak intensity 4 x 10^^ W/cm^, wavelength 1064 
nm, duration T — 10 cycles, and CEP <^ = 0). In panel 



assigned at the initial Mme {t = to = 0) is conserved (^)^ gray-scale color of each trajectory indicates log^o ^ 



over the time evolutio 

We note here that the quantum trajectories can ac- 
tually be propagated without pre-calculating the wave- 
function on the entire simulation space^i^"— . This 
'synthetic' propagation technique has recently attracted 
much attention in view of the development of efficient 
computational schemes to simulate quantum mechanical 
systems with a large number of degrees of freedom. How- 
ever, in the present work, for the sake of focusing on the 
comparison of the Bohmian trajectories with the results 
for the electron probability density we propagated both 
the solution of the TDSE as well as the trajectories. 



III. VISUALIZATION OF THE INTERNAL ELECTRON 
DYNAMICS 

In this section we compare the results from the 
Bohmian trajectory calculations with electron probabil- 
ity densities obtained by integrating the TDSE for the 
ID model interacting with a linearly polarized in- 
tense laser pulse. For this exemplary comparison we have 
chosen the distance between the two protons as i? = 7 
and considered a laser pulse with a peak intensity of 
4 X 10^"^ W/cm^, a wavelength of 1064 nm, a full du- 
ration of T = 10 cycles, and a carrier-to-envelope phase 
of iy9 = 0. We set the soft-core parameter as a = 2.0 
so that the energies of the ground and first-excited elec- 
tronic states of the ID model (—0.519 and —0.491, in- 
cluding the 1/R nuclear repulsion) best reproduce the 
exact values^l (—0.506 and —0.496) for the actual in 
3D space at i? = 7. The initial positions of 100 quan- 
tum trajectories were distributed over — 10 < a; < 10 at a 
regular interval 5z = 0.202, and their weights {wj} were 
determined via Eq. 

In Fig. [T]the Bohmian trajectories (Fig.[TJa)) are pre- 
sented along with the electron probability density (Fig.[T] 



In panel (b), the color code indicates logj^^ where ^ is 
obtained by solving the TDSE (0. 



(b)) as a function of time. For the sake of comparison, 
the electric field of the laser pulse is shown in Fig. [ijc) as 
well. Subject to the intense electric field, the trajectories 
leave the core region (at z w 0) of the molecular ion in 
alternating directions (z — >■ ±oo) at every half cycle of 
the laser field. It is clearly seen that the number of ion- 
izing trajectories increases as the field strength increases 
during the laser pulse. The trajectories liberated from 
the core region show wiggling motions forced by the al- 
ternating electric field of the laser. Due to this quiver 
motion, some of the trajectories, depending on the time 
instants of their release, are driven back to the core region 
and scattered off the protons. Please note that the re- 
sults for the Bohmian trajectories visually agree very well 
with those for the electron probability density: regions of 
large probability density correspond to a large density of 
the Bohmian trajectories. Thus, the Bohmian trajecto- 
ries provide a complete overall picture of the ionization 
process, including the quiver motion and the rescatter- 
ing of the electron in the laser field. This agrees with 
the findings of earlier studies using Bohmian trajectories 

to describe the interaction of atoms with intense laser 
pulses26>28-30^ 

In Fig. [2] we show a detailed view of the electron dy- 
namics in and close to the core region (the protons are 
located at z = ±i?/2 = ±3.5 a.u.) over the central 
field cycle (T/2 - 2n/uj < t < T/2 + 27r/a;) of the laser 
pulse using the electron probability density [Figs, ^a.)- 
mjc)] and the Bohmian trajectories [Figs. [5l^e)-[21[g)]. The 
laser electric field in the same time window is also shown 
in Figs, m^d) andl^^h). In the simulations for Figs. [2l^a) 
and[2{e), the wavefunction masks of cos^/^-shape were 
placed over 270 < |z| < 300 only to avoid the reflection 
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Time (Laser Cycle) Time (Laser Cycle) 

FIG. 2. Close-up view on the time evolution of the electron 
probability density (a-c) and the quantum trajectories (e-g) 
over the single laser cycle at the pulse peak. The laser electric 
field is shown in panels (d) and (h) for comparison with the 
flow of the electron probability density and the quantum tra- 
jectories. In panels (a) and (d), no wavepacket absorber was 
set close to the core region. In panels (b) and (f), an absorber 
was set on one side of the core region at —11 < z < —7. In 
panels (c) and (g), two absorbers were set on both sides of 
the core region at 7 < \z\ < 11. 



of the ionized wavepackets at the boundaries of the large 
simulation box, leaving the electron dynamics in the re- 
gion of z shown in these panels unaffected. Both of the 
electron probability density and the quantum trajectories 
show an ultrafast oscillatory motion around 5 < |z| < 10 
as indicated by, for example, the peaks A-D in Fig.[2je). 
This oscillation involves the ionization and rescattering 
of the electron. To see this we suppressed the rescatter- 
ing effect on the negative z-side by setting a wavefunction 
mask at— ll<z<— 7 for the simulations in Figs. [2jb) 
and[2l^f)i^'^"-2^. As a result, the peak A in Fig.[2]^e) disap- 
pears in Fig. [2{f), indicating that this peak was created 
by the wavepacket driven back from z < and passing 
through the core region toward z > 0. In Figs. [2ljc) and 
[2ljg), the wavefunction masks were set close to the core 
on both sides (at 7 < \z\ < 11) and all the rescattering 
wavepackets were absorbed after the initial ionization. 
The two peaks B and C in Fig. [SJe) are still present in 
Fig. Htg) while peak D disappears, indicating that there 
are actually two bursts of ionization under the present 
conditions. These results confirm the previously reported 
phenomenon of multiple ionization bursts and attosecond 
electron localization^. 

As in Fig. [Tl the quantum trajectories reflect well the 
characteristics of the flow of the electron density in each 
pair of panels in Fig. [2] In addition, the quantum tra- 
jectories in Fig. m^g) clearly elucidates that the electron 



probability moves back and forth between the two pro- 
tons on the ultrafast time scale shorter than a half-cycle 
of the laser field. This is not as obvious in the plot of the 
electron probability density in Fig.[21[c). Thus, the quan- 
tum trajectories are shown to provide further insights in 
the flow of the probability density. 

Please note that the electronic motion in Fig. EJg) be- 
tween the two potential wells created by the protons does 
not necessarily follow the laser-electron interaction po- 
tential Vl. For example, at t = 5 cycles, the oscillating 
laser electric field E{t) is peaked in the negative direc- 
tion (cf. Fig. EJli)), and hence the slope of Vl push- 
ing the electron toward z — > -|-oo becomes maximum. 
Nevertheless, we observe some bound trajectories prop- 
agated in the opposite direction from the proton located 
at z = R/2 = 3.5 to the other one at z = —R/2 = —3.5 
by climbing up the potential Vl- This counter- intuitive 
(and classically forbidden) motion of the electron was no- 
ticed first in the context of coherent control of electron 
localization in dissociating molecule by the Wigner 
representation^. The present results confirm this motion 
and visualize it using Bohmian trajectories. 



IV. ANALYSIS OF THE ELECTRON DYNAMICS USING 
BOHMIAN TRAJECTORIES 

We have seen so far that the Bohmian trajectories 
clearly visualize the transient electron localization and 
multiple ionization bursts within a half-cycle of the laser 
field oscillation. In this section we will now investigate 
the origin of the counter-intuitive motion of the Bohmian 
trajectories. To this end, we focus on the results for 
the intra-molecular electron transfer from one proton to 
the other, obtained by using wavepacket absorbers over 
7 < |z| < 11 in order to eliminate the effect of rescatter- 
ing wavepackets (cf. Figs.[2l^c) and[2l^g)). 



A. Velocity field for the Bohmian trajectories 

As pointed out before, we observe that some of the 
trajectories turn their direction toward z < near the 
peak of the electric field at t = 5 cycles while the classi- 
cal force due to the laser electric field pushes the elec- 
tron in the positive z direction. We now investigate 
this counter-intuitive motion in terms of the velocity 
field, v{x, t) — dS{x, t)/dx, for the Bohmian trajectories. 
Fig.EI^a) shows the phase S{x, t) of the wavefunction and 
Fig. [3fb) the velocity field w(a;,<:) both calculated using 
the same parameters of the field as before. In both Fig- 
ures we superposed the Bohmian trajectories for further 
visualization. 

From Fig.[3]Ja), we can notice that at a given time in- 
stant the phase of the wavefunction is almost constant 
within each of the potential wells. However, the phase 
propagates at different speeds in the two wells. This 
causes a phase gradient around z = and, consequently. 
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B. Phase difference in two-state model 




ill lL S- 

4.75 5 5.25 

Time (Laser Cycle) 

FIG. 3. Analysis of the intra-molecular motion of the quan- 
tum trajectories in terms of the phase of the wavefunction 
S{x,t) and the velocity field v{x,t). (a) Phase S{x,t) and 
the quantum trajectories {zj{t)}. (b) Velocity field v{x,t) 
and the quantum trajectories {zj(t)}. (c) Phase difference 
between the two wells calculated from the ID model (aLR(t), 
blue solid line) and from the approximate solution of the two- 
state model {a^-^^\t), red dashed line), (d) Electric field of 
the laser pulse (wavelength 1064 nm, peak intensity 4 x 10^"^ 
W/cm^, duration T = 10 cycles, and CEP = 0). 



a large absolute value of the velocity field in the region be- 
tween the protons. Over the central half (4.75 < t < 5.25 
laser cycles) of the field cycle shown in Fig.[3l the velocity 
field changes its sign at t = 4.89, 5.0, and 5.10 laser cy- 
cles and, hence, forces the Bohmian trajectories to turn 
their direction at these instants. The change of sign of 
the velocity field corresponds to those instants at which 
the accumulated phase difference between the two wells 
equals a multiple of tt. This interpretation is confirmed 
by the results for the relative phase of the wavefunction, 
defined as 



aLR(0 = arg 



R 
2"' 



arg 



R 



(14) 

and shown as solid line in Fig. ^c). One can clearly see 
that every jump in the relative phase from — tt to 
TT corresponds to an abrupt sign change of v{x,t). We 
may note that this divergence of the v{x,t) is accompa- 
nied by a node in the wavefunction at the same position 
[cf. Fig.HJc)], and that therefore the flux v{x, t)\^{x, t)p 
stays finite. The phase jump between — tt and tt can be 
interpreted as following: If the local phase at one well 
gets more and more retarded (or advanced) from that at 
the other well and if the absolute difference passes tt, then 
the former should be now regarded as more advanced (re- 
tarded) than the latter. 



Next, we show that the phase difference Q!lr can be 
approximated by a simple expression based on a two- 
state model. To this end, we analyze the relative phase 
between the two potential wells in terms of the following 
localized states^"—, 



|L) 
|R) 



1 

'72 
1 

'72 



(15) 
(16) 



where |g) and |u) are the ground and first-excited elec- 
tronic states, respectively, of the ID fixed-nuclei model. 
Without loss of generality, we set the phases of |g) and 
|u) such that |L) and |R) are localized at z = —R/2 and 
z = R/2, respectively. By approximating the state of 
the system in the basis of these two localized states as 
|*(i)) = CL(t)|L) -I- CR(t)|R), the time-evolution of ci,{t) 
and c^ {t) is given by 



CL 

^ dt VcR 



H, 



(21c 



(21o 



with the field-free Hamiltonian 



(21cv) 



Ao /O 1 
2 U 



and the interaction potential 



(21o 



1 

-1 



(17) 



(18) 



(19) 



where Aq is the absolute value of the difference between 
the field-free energies of |u) and |g), and dgu = — (g|z|u) > 
is the transition dipole moment between |g) and |u). 

If the laser-molecule coupling \dguE{t)\ is sufficiently 
strong and/or the laser frequency uj is sufficiently larger 
than the tunnel splitting Ao, we may approximate the 
solution to the two-state TDSE dTT]) by taking y^l^'""') as 
the zeroth order reference Hamiltonian and omitting the 
Ho term2£"— . Such a zeroth order solution can be easily 
obtained as^ 



c^f\t) =exp{-zdg„[A(t) - A{to)]}c^°\to), 
4°)(i) ^cxp {id,^[A{t) - A{to)]}4\to), 



(20) 
(21) 



where the relation of the electric field and vector po- 
tential, Eq. (jH), was used. Then, the phase difference 
between the two wells can be approximated by 



(21ev) / ,x 

"lr W = arg 



,(0) 



arg 



,(0) 



(22) 



By substituting Eqs. ^ and dH]) into Eq. dH]), and 
noting that A(to) = at the beginning of the laser pulse 
at i = to = 0, we obtain 



a 



(21cv) 
LR 



(t) = -2dguAit) + 



a 



(21cv) 
LR 



(to). 



(23) 
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The phase differenee calculated by this expression is plot- 
ted in Fig. ^c) as the red dashed line. We can see that 
the two-state model closely reproduces the exact value 
(blue solid line) calculated numerically by Eq. p^ . 

The condition, a^^^\tturn) = nir, n G Z, for the time 
instant fturn at which the velocity field between the two 
protons changes its sign is then given by 



Aitt 



(21cv) 



(to) 



2d, 



(24) 



This expression has a similar form as the condition for 
the time instant tioc of maximum electron localization, 



A(tloc) = 



X 



'2(1 a; Yi 



(25) 



with m li derived previously^ based on the series ex- 
pansion of the Floquet states for the two-state model^. 
In fact, these two expressions are identical at the limit of 
long laser pulse duration where the mixing angle x of the 
two Floquet states reduces to the initial phase difference 



a 



(21ev) 
LR 



(to) 



C. Origin of the counter-intuitive electron motion 



The result of the two-state analysis in the previous 



subsection, in which we consider Vj 



(21ev) 



as the zeroth 



order Hamiltonian while neglecting the tunnel hopping 
term Hq, provides us with an intuitive picture. Please 
note that the transition dipole moment has the asymp- 
totic form dgu ~ i?/2 at large F^, and the difference of 
the diagonal elements of l/^^icv) hence, approximately 
equal to RE{t), which is the difference of the electric po- 
tential induced by the laser field between the two wells. 
Using this asymptotic form of dgu, the phase difference 
in Eq. can be rewritten as 



(21ev) 



it) 



dt' RE{t') + Jl^"'\to) 



to 



(26) 



This expression elucidates that the origin of the phase 
difference between the two potential wells, and hence the 
velocity field, is the difference of the electric potential 
energies between the two wells induced by the laser light. 

This interpretation of the electronic dynamics based 
on the zeroth order two-state analysis predicts that the 
motion of the electron should be still counter-intuitive 
even at a relatively low intensity at which only one lo- 
calization per half-cycle is predicted by Eq. ([M|) or (PSI) . 
This is demonstrated in Fig. |4] in which the ID model 
was initially prepared in the ground state, and a laser 
pulse of peak intensity 4 x 10"'^^ W/cm^ was applied. The 
other parameters were the same as above. The initial 
phase difference C(^-i^^\to) is zero due to the choice of the 
initial state as |g), and the velocity field around z = ev- 
olcs in time as —A{t) according to the two-state analysis 




4.75 5 5.25 

Time (Laser Cycle) 

FIG. 4. Time evolutions of (a) the phase of the wavefunction 
and (b) the velocity field are plotted along with the quan- 
tum trajectories. The system was initially prepared in the 
ground state |g). Panel (c) shows the phase difference be- 
tween the two local potential wells calculated from the TDSE 
solution (blue solid line) and from the two-state model (red 
dashed line). Panel (d) shows the electric field of the laser 
pulse (wavelength 1064 nm, peak intensity 4 x lO^'^ W/cm^, 
duration T = 10 cycles, and CEP (p ^ 0). 



(Fig-Sic)). Due to the phase-lag between A{t) and E{t), 
the quantum trajectories are navigated from z = R/2 to 
z = —R/2 during the 4.5 < t < 5, for example, while 
the electric force —E{t) points in the opposite direction 
during 4.75 < t < 5. As a consequence, the trajectories 
are always accumulated at the upper potential well, in 
contradiction to our classical intuition. 

If, instead, the system is initially prepared in the first 
excited state |u), the initial phase difference a'^^^\to) = 
— TT. Due to this offset, the quantum trajectories (and 
the electron probability density) are navigated toward 
the lower potential well in this case. In fact, we can see 
that this is the case for the results presented in Fig. [5] 
where we applied the same laser pulse as used for the 
results shown in Fig. |4] but prepared the system in the 
ID model of in the first excited state. 

We may finally ask if there is a parameter regime in 
which our intuitive (classical) picture that the electron 
should move in the direction of the electric force is always 
recovered. As we mentioned above, our two-state analy- 
sis is based on the assumption that y^^^^'"^) is the domi- 
nant term, i.e. that the laser-molecule coupling is strong 
and/or the photon energy is sufficiently larger than the 
tunnel splitting. Please note that the photon energy at 
the wavelength of 1064 nm is a; = 0.043, whereas the tun- 
nel splitting is Aq = 0.028 for the ID model at i? = 7. 
If we decrease the photon energy as well as the laser in- 
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FIG. 5. Same as Fig. |4] except that here the system was 
initially prepared in the first excited state |u). 
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FIG. 6. Same as Fig. |4]except that here a laser pulse with long 
wavelength (5000 nm) and low intensity (4 x 10^° W/cm^) 
(duration T = 10 eye, CEP — 0) was applied. 



tensity, the two-state analysis is no longer applicable. In 
fact, the results in Fig. [6] show that at a wavelength of 
5000 nm and an intensity of 4x 10^° W/cm^ our two-state 

analysis breaks down, since the phase difference cx^^^^ (t) 
(red dashed line) predicted by the two-state analysis de- 
viates qualitatively from the value of aLR(i) (blue solid 
line) calculated from the TDSE solution (Fig. [SJc)). We 
can see that in this parameter regime, in which the con- 
ventional weak field perturbation theory may be applied. 



the motion of the quantum trajectories (and the electron 
probability density) follows the electric force of the laser 
field, and the intuitive (classical) picture is indeed recov- 
ered. 



V. CONCLUSIONS 

We have presented an analysis of the attosecond elec- 
tron dynamics in hydrogen molecular ion driven by an 
intense laser pulse in terms of Bohmian trajectory calcu- 
lations using a ID model. Results of these calculations 
are found to be in good overall agreement with those of 
numerical simulations of the TDSE. Recently predicted 
phenomena such attosecond transient electron localiza- 
tion and multiple bursts of ionization within a half cycle 
of the laser pulse are clearly represented by the Bohmian 
trajectories. Further analysis let us identify the origin of 
the sometimes counter-intuitive motion of the Bohmian 
trajectories as due to the time evolving phase difference 
of the wavefunction between the two potential wells in- 
duced by the electric potential of the laser field. We were 
able to predict the time instants at which the trajectories 
change their directions within a simple two-state model. 
Following our analysis we showed that exposed to an in- 
tense laser field the electron dynamics in the hydrogen 
molecular ion often does not follow the (classical) force 
of the laser electric field. Our classical expectations are 
however recovered in the perturbative weak-field limit in 
which the photon energy and intensity of the laser field 
are both sufficiently small. 
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